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The microstructure of coagulated colloidal particles, for which the inter-particle potential is de- 
scribed by the DLVO theory, is strongly influenced by the particles' surface potential. Depending on 
its value, the resulting microstructures are either more "homogeneous" or more "heterogeneous", at 
equal volume fractions. An adequate quantification of a structure's degree of heterogeneity (DOH) 
however does not yet exist. In this work, methods to quantify and thus classify the DOH of mi- 
crostructures are investigated and compared. Three methods are evaluated using particle packings 
generated by Brownian dynamics simulations: (1) the pore size distribution, (2) the density fluc- 
tuation method and (3) the Voronoi volume distribution. Each method provides a scalar measure, 
either via a parameter in a fit function or an integral, which correlates with the heterogeneity of 
the microstructure and which thus allows for the first time to quantitatively capture the DOH of 
a granular material. An analysis of the differences in the density fluctuations between two struc- 
tures additionally allows for a detailed determination of the length scale on which differences in 
heterogeneity are most pronounced. 

PACS numbers: 81.05.Rm, 61.43.-j, 82.70.Gg 



INTRODUCTION 



Colloidal particle packings are suitable model systems 
for the study of the structural properties of granular ma- 
terials below the random loose packing limit. For such 
systems the gravitational force is negligible in compari- 
son to the van der Waals force, electrostatic repulsion or 
Brownian motion P, In the present study, we par- 
ticularly focus on systems, for which the local arrange- 
ment of the particles is the only variable, as opposed 
to variations in the volume fraction or the particle size 
distribution for example. Commonly, these microstruc- 
tures are referred to as "more homogeneous" or "more 
heterogeneous" , which either designates a structure pre- 
senting a rather uniform distribution of the particle po- 
sitions or one having locally denser regions and therefore 
larger voids. These qualitative terms may be intuitive, 
however, they do not allow for a sound scientific quantifi- 
cation of the structure's degree of heterogeneity (DOH), 
which does not yet exist. In this paper, three methods 
providing the means for such a quantification are pre- 
sented, analyzed and compared. These methods permit 
for the first time to explicitly capture the DOH of a par- 
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tide packing in form of a quantitative, scalar measure. 

Experimentally, the reproducible generation of col- 
loidal particle packings possessing a specific DOH is 
achieved by the use of an in-situ enzyme-catalyzed 
destabilization method (Direct Coagulation Casting 
(DCC), 0,111). For volume fractions between 0.2 and 0.6, 
DCC allows for the coagulation of electrostatically stabi- 
lized colloidal suspensions to stiff particle structures by 
either shifting the pH of the suspension to the particles' 
isoelectric point or by increasing the ionic strength of the 
suspension without disturbing the particle system. Shift- 
ing the pH leads to " more homogeneous" microstructure 
through diffusion limited aggregation while increasing 
the ionic strength results in "more heterogeneous" mi- 
crostructures via reaction rate limited aggregation. The 
differences in heterogeneity have been observed using var- 
ious techniques such as diffusing wave spectroscopy 0], 
static light transmission Q or cryogenic scanning elec- 
tron microscopy ^Sj. 

Computationally, microstructures with different DOH 
were successfully reproduced using Brownian dynam- 
ics simulations (BD) [1, 0, [!]■ Slices of three particle 
layer thickness through a "homogeneous" and a "het- 
erogeneous" BD-microstructure of identical volume frac- 
tion nicely demonstrate the differences between the mi- 
crostructures (Fig. [T]). The microstructure on the left 
presents a rather uniform distribution of the particle po- 
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FIG. 1: Slices through a "homogeneous" (left) and a "hetero- 
geneous" (right) particle structure with the same volume frac- 
tion of 0.4 resulting from Brownian dynamics simulations 01 
(shea thickness: three particle diameters; particle diameter: 
0.5 fj,m). 



sitions over the whole slice while the microstructure on 
the right presents locally more densely packed particles 
and consequently larger voids. Both structures have an 
identical overall volume fraction of 0.4. 

In preceding works, various characterization methods, 
such as the radial pair correlation function [7] , the bond 
angle distribution function JTj, the triangle distribution 
function 0] and the Minkowski functionals in conjunction 
with the parallel-body technique ^ were applied to sets 
of microstructures generated by BD simulations Q . 

The pair correlation function quantifies the amount of 
structural rearrangement during the coagulation. Its use- 
fulness regarding a distinction between structures with 
different DOH however is rather limited as the differences 
between peaks corresponding to characteristic particle 
separation distances are relatively small 7]. The main 
advantage of the pair correlation function is its exper- 
imental accessibility through scattering techniques such 
as SESANS 0. 

The bond angle distribution function and the triangle 
distribution function give further information on the lo- 
cal building blocks of the particle network Q . Particular 
features, as for example peaks in the respective distri- 
bution function allow distinguishing between "more ho- 
mogeneous" and "more heterogeneous" microstructures. 
However, as in the case of the pair correlation function, 
the differences between structures with different DOH are 
small for both, the bond angle and the triangle distribu- 
tion function. 

The analysis using the Minkowski functionals in con- 
junction with the parallel body technique supplies addi- 
tional information on the structure's morphology resolv- 
ing microstructural differences on a length scale limited 
by the largest pore size @. This method is computa- 
tionally intensive and the extension to arbitrary particle 
shapes is difficult. 

Gearing towards a possible correlation between mi- 
crostructure and mechanical properties "homogeneous" 
and "heterogeneous" microstructures have recently been 
analyzed in terms of load bearing sub-structures: Firstly, 



regions of closely packed particles and secondly, quasi- 
linear chains of contacting particles [lH. The locally 
closed packed regions were analyzed using the common 
neighbor distribution in conjunction with the dihedral 
angle distribution. Both methods only showed minor 
differences between a "homogeneous" and a "heteroge- 
neous" microstructure. In particular, practically the 
same number of triangles and regular tetrahedrons were 
found in both structures. Quasi-linear arrangements of 
contacting particles were quantified using the straight 
path method, revealing significant microstructural dif- 
ferences between "homogeneous" and "heterogeneous": 
approximately twice as many paths of length longer or 
equal to four particles and three times as many paths of 
length longer or equal to five particles where found in the 
" heterogeneous" microstructure. 

Despite the multitude of microstructural characteri- 
zation methods that have been applied to colloidal mi- 
crostructures possessing various DOH, a useful quantifi- 
cation of the microstructures' heterogeneity in form of a 
scalar measure is lacking. In the present study statisti- 
cal measures allowing for a quantification of a structure's 
heterogeneity are provided. The following three methods 
aiming at this quantification are discussed. 

Firstly, the exclusion probability [l^l estimates the 
pore size distribution by randomly probing the pore 
space. In [13[ this method was applied to very dilute 
simulated colloidal systems and allowed for a clear dis- 
tinction between particle gel networks with varying tex- 
tures. The same method, termed spherical contact dis- 
tribution function was used to investigate the pore 
size distribution of dense sphere packings as a function of 
the particle size distribution and the packing generation 
algorithm. 

Secondly, the density fluctuation method considers the 
fixed particle centers as a point process, and it statisti- 
cally analyses the fluctuations of the particle center den- 
sity as a function of length scale. Used comparatively, 
this method further allows for a detailed analysis of the 
length scale, on which two structures present the largest 
differences in terms of heterogeneity. 

Finally, the Voronoi volume distribution [l^ is used 
to quantify the distribution of the free volume of our 
particle packings. In ^6|] the packing of cohesive par- 
ticles resulting from simulations using the discrete ele- 
ment method [17] with volume fractions between approx- 
imately 0.2 and 0.6 was analyzed. The distribution of 
Voronoi volumes was shown to broaden with decreasing 
volume fraction. In [l8| the Voronoi volume distribution 
was determined for a large set of experimental and nu- 
merical data covering a wide range of volume fractions. 
The various distributions were shown to follow a so-called 
/c-gamma function, which was deduced by means of a 
statistical mechanics approach. The parameter k charac- 
terizing the shape of the curve was found to depend very 
sensitively on the structural organization of the particles. 

To the authors' knowledge, none of these methods have 
yet been applied to particle structures, for which the 
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TABLE I: DLVO Potential Parameters. 



Parameter 


Symbol Value 


Hamaker constant of AI2O3 in H2O 


Ah 


4.76 X 10"^" 


Particle diameter 


do 


5 X 10"'^ m 


Relative dielectric constant of H2O 


€r 


81 


Surface potential 


*0 


- 15 mV 


Absolute temperature 


T 


293 K 


Valency of ions 


z 


1 


Inverse Debye screening length 




10** m"^ 



DOH is the only variable, in opposition to a varying vol- 
ume fraction or particle size distribution for example. 



II. MATERIALS AND METHODS 

In the following, the structure characterization meth- 
ods employed in this work are presented: the pore size 
distribution, the density fluctuation method and the dis- 
tribution of Voronoi volumes. These methods are eval- 
uated in terms of their ability to quantify the DOH 
of microstructures generated by previous BD simula- 
tions 0, In these simulations, the DLVO-theory [l^ 
was used to describe the inter-particle potential V^^^" 
given by the sum of an attractive van der Waals term 
yvdw j|gq^ [jjj g^j-^j electrostatic repulsion term V^' 
(Eq.O.Thus, y^'™ = V"'^'" + with. 
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FIG. 2: Particle structure (gray particles) with a test particle 
inserted at position P. 



consist of 8000 particles and are contained in a simula- 
tion box with periodic boundary conditions. In particu- 
lar, the particle interpenetration is much smaller than the 
length scale of the heterogeneities analyzed in this study. 
Less than 0.1% of all contacts present an interpenetration 
above 1.0% of the particle diameter do with a maximum 
interpenetration of 1.3% c?o. Please refer to 0, Q for a 
more detailed description of the BD-simulations. 
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respectively. The DLVO parameters are summarized 
in Table n 

The heterogeneity of the final microstructure was 
shown to be closely related to the presence and depth of 
the secondary minimum in the inter-particle potential, 
which, for fixed values of the Debye screening length, 
essentially depends on the particles' surface potential 

The microstructures analyzed in this work are labelled 
according to the surface potential VPo used during their 
generation. In particular, the following 5'o-values are 
used: mV, 12 mV, 13 mV, 14 mV and 15 mV. Addi- 
tionally, these final microstructures are compared to the 
initial microstructure, representing a stabilized suspen- 
sion in which the inter-particle potential is purely repul- 
sive. The volume fraction is fixed at 0.4, the monosized 
particles have a radius ro = 0.25 /zm, all microstructures 



A. Pore Size Distribution 

The pore size distribution is estimated following the 
approach described by Torquato et al. [l2| using the ex- 
clusion probability Ev{r). Ev{r) is defined as the prob- 
ability of inserting a "test" particle of radius r at some 
arbitrary position in the pore space of a system of N par- 
ticles. This is schematically represented in Fig. [5] using a 
set of particles of radius tq (gray) with a test particle of 
radius r inserted at position P. 

In order to estimate Ey {r) , a statistically large number 
of points is randomly placed in the pore space of a given 
microstructure and the distance to the closest particle 
surface is determined. 

In P the relation between Evij) and the Minkowski 
functional Wi{r) was described. The Minkowski func- 
tional in conjunction with the parallel body technique 
considers the point process given by the fixed particle 
centers. Generally, in three dimensions there are four 
functionals Wi, where i = 0, with d the spatial di- 
mension, corresponding to the volume, surface, aver- 
age mean curvatures and connectivity. In particular. 
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Wi{r') = ^ /^^(r') ^'^ calculates the surface of the union 
of spheres located at the particle centers in dependence 
of their radius r' — rQ+r. Schematically, Wi{r') is shown 
in Fig. [21 Wi{r')dr' is the volume between the distance 
r' and r' + dr' . The probability of placing an uniformly 
drawn random test point at a distance re [r' ,r' + dr'[ is 
proportional to the volume delimited by r' and r' + dr' 
and therefore P{r' < ro + r < r' + dr') = Ev{r)dr' cx 
Wi{r')dr' , which results in Ev{ro + r) oc Wi{r'). Thus, 
given a statistically large number of test points Ev{r) 
provides a means to estimate Wi(r'), which has the ad- 
vantage of being computationally less intensive. A simi- 
lar Monte Carlo integration is usually performed to cal- 
culate Wo(r'). 



B. Density Fluctuations 

The density fluctuation method quantifies the local 
fluctuation of particle centers by subdividing the struc- 
ture into smaller parts and measuring the average value 
and the standard deviation of the particle center density. 
Practically, this is done by dividing the structure into 
cells using a cubic grid, where ric is the number of cells 
along one dimension with ric — 2, .., njy*"^, n™"^ being 
the maximum number of cells under consideration (along 
one dimension). The density fluctuation method consists 
in determining the average number of particle centers per 
cell Eppc and its standard deviation Uppc as a function of 
ric and then calculating the relative fluctuations . 

Additionally, the density fluctuation method is applied 
comparatively allowing for a determination of the length 
scale on which two structures show the largest differences. 
Therefore, the difference A(nc) given in Eq. ([3]) is calcu- 
lated between two structures i and j . 



' ppc 



El 



E^ 



(ric) 
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If both structures have the same number of particles 
and identical volume fractions, which is the case for the 
microstructures investigated in this work, then Epp^ = 
E'ppc = Eppc and Eq. ^ yields Eq. (U). 



is a polyhedron whose interior consists of all points in 
space that are closer to the center of particle i than to 
any other particle center [l^. The Voronoi tessellation 
thus divides the volume containing a set of particles into 
a set of space-filling, non-overlapping and convex polyhe- 
drons. In this work, the Quickhull algorithm [2l[ is used 
to compute the volumes of the Voronoi polyhedrons. The 
distribution of the Voronoi volumes describes the devia- 
tion of a structure from a perfect crystalline packing, in 
which case all particles occupy the same volume and the 
Voronoi volume distribution thus is a delta function. The 
minimum volume of a Voronoi cell Vmin is achieved for 
a regular close packing with Vmin — l-325T4pft,ere, where 
Vsphere is the volumc of a particle. The difference be- 
tween a particle's Voronoi volume Vi and Vmin is termed 
the Voronoi free volume V/ — Vi — Vmin- 

The distribution of the Voronoi free volume was found 
to follow gamma-distributions: Kumar and Kumaran for 
example have shown that the free volume distribution 
of hard-disk and hard-sphere systems are well described 
using a two-parameter and a three-parameter gamma- 
distributions (2^ . 

Aste et al. have deduced the two-parameter gamma- 
distribution using a statistical mechanics approach p^ . 
The so-called fc-gamma distribution given by Eq. ([5]) was 
found to agree very well with a large number of exper- 
iments and computer simulations over a wide range of 
packing fractions. 



/(W,fc) 



{Vf) 



k-l 



r(fc) 



■exp( 



Vf 



(5) 



The mean Voronoi free volume V^ is a scaling param- 
eter given by y - Vmin = (1/$ - 1.325) Vsphere, where $ 
is the volume fraction. The free parameter k characteriz- 
ing the shape of the curve depends very sensitively on the 
structural organization of the particles and corresponds 
to the specific heat in classical thermodynamics. Empir- 
ically, k can be computed using k — -^^3^, where ay is 

the variance of the free volume distribution. In particu- 
lar, Eq. ([5]) was shown to hold for systems at statistical 
equilibrium as well as for systems out of equilibrium . 

III. RESULTS AND DISCUSSION 



A(nc) 



En 



(ric) 
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A. Pore Size Distribution 



A is plotted against the cell's edge length lc{nc) = 
LBox/nc normalized by the particle diameter do, where 
Lbox is the side length of the cubic simulation box. 



C. Voronoi Volume Distribution 

Formally, for a set of monodispersed spherical parti- 
cles, the Voronoi volume Vi associated with a particle i 



The pore size distribution of the various microstruc- 
tures is depicted in Fig. [3] in terms of the exclusion prob- 
ability Ev{rp) as a function of the pore radius rp nor- 
malized by the particle radius rp. A set of 10^ random 
test points [11] placed in the structures' pore space was 
used to estimate Ev{rp). 

All curves for the final microstructures (^'o = - 
15 mV) decrease monotonically towards increasing pore 
sizes indicating a decreasing probability of finding larger 
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FIG. 3: Pore size distribution for the initial microstructure 
(stable suspension) and final microstructures (^'o = 0-15 
mV). B-spline curves serve as guide to the eye. 



FIG. 4: Probability of finding pores with a radius rp larger 
than r vs. r normalized by the particle radius ro for the vari- 
ous microstructures (symbols) and corresponding fits using a 
complementary error function (solid lines). 



pores. A particular behavior is found for the stable sus- 
pension where the exclusion probability increases with 
increasing pore size up to a pore radius of O.lSrQ. This 
is due to the repulsive inter-particle potential in the case 
of the stable suspension where consequently, the particles 
are not in contact. The pore diameter, at which the max- 
imum in the exclusion probability is found, corresponds 
to the average surface-to-surface distance between neigh- 
boring particles. 

Remarkably, the various curves for the final mi- 
crostructures in Fig. [3] intersect at approximately the 
same point defining a characteristic pore size r^, found 
at 0.65ro. The probability of finding a pore with a radius 
below Tp decreases for increasing values of while pores 
with a radius above are found with higher probability 
towards increasing ^^q. Indeed, the probability of finding 
pore radii larger than approximately l.lrg is negligible in 
the case of the " most homogeneous" microstructure with 
^'o = mV while pore radii up to 2.4ro are found in the 
"most heterogeneous" microstructure with — 15 mV. 

Using the results obtained for the exclusion probability 
Evifp), the cumulative probability P{rp > r) of finding 
pore radii larger than r was calculated using P(rp > r) = 
J2r'>rp Ev{r'). The results are shown in Fig. [D 

P{rp > r) decreases monotonically for all microstruc- 
tures. The fastest decrease is found for the stable sus- 
pension. With increasing '^q the decrease of P{rp > r) 
is slower. Comparing the "most and the least hetero- 
geneous" microstructure, with ^Pq = 15 mV and ^'o — 
mV, respectively, there is a 1.7 times higher probabil- 
ity of finding pores larger than 0.5ro. Towards larger 
pore radii, the probability ratio increases: finding pores 
with a radius larger than 0.75ro and l.Oro is 5.1 times 
and 60 times, respectively, more probable in the "het- 
erogeneous" than in the "homogeneous" microstructure. 
Fig. 0] further shows the fit of P(rp > r) using a comple- 
mentary error function given by Eq. 



P{rp > r) = 1 - erf 



r/ro 



i\/2 



(6) 



The error function is defined as the cumulative Gaus- 
sian distribution: erf(a;) = JJ^ exp (— z'^) dz. Pa- 
rameter a is the standard deviation, i.e. the width of the 
corresponding Gaussian distribution and h is the location 
of its maximum, i.e. the most probable pore to particle 
radius ratio. 

Table HIl summarizes the fit parameter a and b obtained 
for the various microstructures analyzed in this work and 
the corresponding correlation coefficients R^, which for 
all fits are very close to 1 and thus indicate a good fit. 
Parameter a is smallest for the initial microstructure and 
increases towards increasing values of ^'o- The increas- 
ing values of a reflect the slower decrease of the curves 
in Fig. [4] and hence the broadening of the distributions 
towards increasing DOH. The values found for b decrease 
with increasing ^E'o representing a shift of the maximum 
in the Gaussian distribution shown in Fig. [31 

In the following expression for Ev{r) was found 
for a statistically homogeneous microstructure of impen- 
etrable spheres: Ev{r) = (1 — $) exp(P3(r, $)), where 
$ is the volume fraction and P3 is a third degree poly- 
nomial function in r. This function can be interpreted 
as a corrected Gaussian distribution which is nicely ap- 
proximated by Eq. ^ as well (R^ — 0.9962), yielding 
a = 0.347 and b = 0.0204. The DOH of this theoretical 
structure thus lies between the stable suspension and the 
most homogeneous final microstructure with '5 = mV. 

An alternative to the quantification of a structure's 
heterogeneity by means of the fit parameter a is the cal- 
culation of the integral over the cumulative pore size dis- 
tribution. This scalar measure has the advantage of being 
statistically more robust. It is also more general in the 
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TABLE II: Fit parameters a and fe, R^-values and integrals 
/ps and Idf for the various microstructures. 



b (10- 



stable suspension 0.2651 2.042 0.9952 0.215 20.29 

*o = mV 0.3752 0.9013 0.9987 0.291 22.01 

*o = 12 mV 0.4053 0.3136 0.9994 0.311 22.31 

*o = 13 mV 0.4127 0.1225 0.9996 0.315 22.44 

*o = 14 mV 0.4736 -1.415 0.9996 0.352 23.05 

*o = 15 mV 0.5377 -3.756 0.9972 0.388 23.71 



sense that it is applicable even when the fit with a com- 
plementary error function does not yield good results. 
The integral over the cumulative pore size distribution is 
labelled Ips and is given by Eq. ([7]): 



Ips = / P{rp >r)—^ — Y\ P{rp > n) 
Jr>o ro ro 



(7) 



The second equality accounts for the discrete case, 
where 5r is the radial resolution of the empirical pore size 
distribution. The /pg-values for the various microstruc- 
tures are summarized in Table |TT1 In our case, in which 
the data can nicely be fitted using Eq. ([6]), the fit pa- 
rameter a is proportional to Ipg: a/Ipg = 1.31 ± 0.05. 
Thus, a quantification of the DOH by means of a or Lps 
is equivalent. 



B. Density Fluctuations 

The density fluctuations for the various microstruc- 
tures are shown in Fig. [51 Over the whole range of 
grid spacings, the fluctuations are smallest for the sta- 
ble suspension and increase for increasing values ^'q. For 
Tie > 34, which corresponds to a grid spacing of 0.64 
particle diameter, the density fluctuations of the various 
microstructures are equal. Idf given in Eq. ^ provides 
an integral measure of the heterogeneity similar to Ips in 
the previous section, however accounting for the discrete 
variable rir. 



n,<34 ^PP" 



(8) 



The Idf values for the various microstructures summa- 
rized in Table Hi] continuously increase towards increas- 
ing values of and thus measure the DOH of the mi- 
crostructures. 

In the following, two sets of comparisons are per- 
formed: Firstly, the various final microstructures are 
compared to the stable suspension. This comparison 
quantifies the length scale on which structural rearrange- 
ments take place during the coagulation. Secondly, the 
final microstructures with > mV are compared to 
the "most homogeneous" microstructure with = 
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FIG. 5: Relative density fluctuations for the various final mi- 
crostructures and the stable suspension as a function of grid 
spacing. 



mV. This set of comparisons reveals the length scale on 
which variations in heterogeneity are most pronounced. 

The comparison of the various final microstructures to 
the initial, stabilized microstructure is shown in Fig. [5] 
in terms of A(rtc) as given by Eq. ^ where superscripts 
i and j correspond to a final microstructure and to the 
initial microstructure, respectively. A(nc) is shown as 
function of the grid spacing / normalized by the particle 
diameter do = 2ro. 

All curves present an identical behavior, which essen- 
tially consists of three successive peaks with decreasing 
height towards a larger grid spacing. The location of the 
first, second and third peak is slightly above one, at two 
and at three particle diameters, respectively. The height 
of the individual peaks increases for increasing values 5*0. 

More precisely speaking, the first peak is found at 
1.09(io for all final microstructures in comparison to the 
stable suspension. This grid spacing corresponds to a cell 
number of 8000 and therefore to the case where the aver- 
age number of particles per cell is exactly one. This case 
is best reproduced for the stable suspension as for the fi- 
nal microstructures the standard deviation is roughly 20 
to 27% larger. Physically, this peak is explained by the 
transition of the inter-particle potential from repulsive 
to attractive. Indeed, the repulsive potential in the case 
of the stable suspension causes all particles to occupy 
approximately the same volume as will be confirmed in 
Sec. nil CI by means of the Voronoi volume distribution. 
The switching of the inter-particle potential from repul- 
sive to attractive causes the particles to form contacts 
resulting in an average particle separation of one par- 
ticle diameter, which is smaller than the grid spacing 
of l.OQdo. This increases the probability of finding cells 
that are cither empty or contain more than one particle 
and thus the standard deviation of the average number 
of particles per cell is increased. 

The peaks at a grid spacing of approximately two and 
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FIG. 6: Relative difference between the density fluctuations 
of tiie various final microstructures and the stable suspension 
as a function of grid spacing. B-spline curves serve as guide 
to the eye. 



FIG. 7: Relative difference between the density fluctuations of 
the various microstructures with ('I'o > mV) and the "most 
homogeneous" microstructure ('I'o — mV) as a function of 
grid spacing. B-spline curves serve as guide to the eye. 



three particle diameters are considerably less pronounced 
than the peak at 1.09do. In particular, the differences 
between the various final microstructures are larger than 
for the first peak. These differences will be elaborated in 
more detail in the following. 

The comparison between the final microstructures with 
^'o > mV and the " most homogeneous" microstructure 
with ^'o = mV is shown in Fig. [T] Here, superscripts i 
and j (Eq. ([¥])) correspond to one of the microstructures 
with > mV and to the microstructure with '^q — 
mV, respectively. Over the whole range of grid spacings, 
the differences between the density fluctuations increase 
towards higher values of ^'o- This behavior correlates 
very well with the increase of porosity for increasing 
as already observed in the previous section. Addition- 
ally, Fig. [7] reveals that the largest differences in terms 
or particle density between the "most and least hetero- 
geneous" microstructure (vPo = 15 mV and ^'o = mV, 
respectively) are found on a length scale between 1.3 and 
2.2 particle diameters. 



Voronoi Volume Distribution 



As stated in Sec. Ill CI the distribution of Voronoi vol- 
umes Piv) describes the deviation of a given structure 
from a perfectly crystalline packing, for which P(i') is a 
delta function and all particles occupy the same volume. 
For random particle structures P(i^) broadens and as will 
be shown in the following, the width of the distribution 
can be interpreted as the heterogeneity of a structure. 

We've calculated P(j^) for the stable suspension and 
the various final microstructures as a function of = 
the free volume normalized by the mean free 
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13 mV 

14 mV 

15 mV 
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volume (Fig. [51 symbols). The distribution found for the 
stable suspension is significantly narrower in comparison 



FIG. 8: Voronoi volume distribution Pif) for the various mi- 
crostructures (symbols) and corresponding fits using the k- 
gamma distribution (lines). 



to the final microstructures. In the case of the final mi- 
crostructures, P{ry) broadens with increasing value of 
indicating that larger fluctuations in Voronoi volumes are 
found with increasing heterogeneity. 

The various microstructures were fitted using the k- 
Gamma distribution given in Eq. ©. The resulting 
curves are shown in Fig. [8] (lines) and the correspond- 
ing k- and R^-values are summarized in Table ITTll k de- 
creases with increasing heterogeneity and can therefore 
be used as a measure for the DOH of the microstruc- 
tures. The R^-values close to one indicate good fits. 
In particular, a very good fit quality was achieved for 
the stable and the final microstructures up to — 14 
mV. The R^-value for the ^'o = 15 mV microstructure 
is lower. Indeed, Fig. [8] shows that for the \l/o — 15 mV 
microstructure the agreement between the measured dis- 
tribution and the fit curve for small Voronoi volumes is 
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TABLE III: fc-Gamma fit results. 



k 

stable suspension 62.3 0.990 

*o = mV 8.6 0.995 

*o = 12 mV 6.5 0.996 

*o = 13 mV 6.0 0.996 

*o = 14 mV 4.0 0.994 

*o = 15 mV 2.8 0.972 



not as good as for the other curves. This might be re- 
lated to the fact that during the generation of the = 15 
mV microstructure the energy barrier between primary 
and secondary minimum in the inter-particle potential 
was largest. This resulted in a few particle contacts still 
trapped in the secondary minimum (roughly 7% of the 
physical contacts). Particles trapped in the secondary 
minimum have an inter-particle distance of 2.16ro instead 
of 2ro upon complete coagulation, which may be a rea- 
son for the reduced fit quality towards smaller Voronoi 
volumes. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have analyzed three distinct mi- 
crostructural characterization methods. Using these 
methods, scalar measures were introduced, which for the 
first time allow quantifying the DOH of particle packings. 

• The exclusion probability gives an estimate of the 
pore size distribution by a random probing of 
the pore space. The "more heterogeneous" mi- 
crostructures present a considerably broader pore 
size distribution with a significantly longer tail than 
the distribution for the "more homogeneous" mi- 
crostructures. In particular, a continuous broad- 
ening is found with increasing heterogeneity. The 
cumulative exclusion probabilities were shown to 
follow error functions with parameter a reflect- 
ing their width and thus measuring the structures' 
DOH. Fit parameter a increases with increasing 
heterogeneity. 

• The density fluctuation method statistically ana- 
lyzes the particle center density in dependence of 
the sampling length scale. The relative density fluc- 
tuation as function of the grid spacing presents a 
clear dependence on the heterogeneity of the mi- 
crostructure: Over the whole range of grid spac- 
ings, the stable suspension exhibits the smallest 
density fluctuations. These fluctuations increase 
towards increasing values of and thus increasing 
DOH, which is nicely reflected by increasing values 
of the integral measure Idf- 

An examination of the differences between the den- 
sity fluctuations of two structures was found to 
be particularly useful as it allows determining the 




FIG. 9: Interdependence between the measures of the degree 
of heterogeneity for the various microstructures: Idf and 1/k 
as a function of a. 

length scale on which the structures present the 
largest differences in heterogeneity. In the case 
of the "most heterogeneous" microstructure the 
largest differences in comparison to the "most ho- 
mogeneous" one are found on a length scale be- 
tween 1.3 and 2.2 particle diameters. 

• The Voronoi volume distribution of the stable sus- 
pension is very narrow in comparison to the final 
microstructures, for which the distribution broad- 
ens with increasing heterogeneity. The various 
Voronoi volume curves were shown to follow k- 
gamma distributions. Parameter fc, reflecting the 
width of the distribution and thus the structure's 
DOH, decreases with increasing heterogeneity. 

The behavior of the three parameters a, k and Idf is 
summarized in Fig. [5] showing Idf (left scale) and 1/k 
(right scale) as a function of a. The solid lines suggest a 
pairwise affine dependence between Idf , 1/k and a. Thus, 
as far as the quantitative characterization of the DOH of 
the microstructures considered in this work is concerned, 
all methods, the pore size distribution, the density fluc- 
tuation method and the Voronoi volume distribution, can 
be considered as equivalent in the sense that the knowl- 
edge of one parameter permits to determine the others. 
However, parameter k reflects changes in the DOH more 
sensitively than a or Idf. Indeed, the normalization of fc, 
a and Idf with respect to their maximum values reveals 
that parameter k covers the interval [0.04, 1.0]. This in- 
terval is significantly larger than the normalized ranges of 
a and Idf, which are [0.5, 1.0] and [0.86, 1.0], respectively. 

The interrelation between the three structural charac- 
terization methods can be understood as follows: The 
probability of placing a random point used for the de- 
termination of the pore size distribution into the free 
Voronoi volume of a particle is proportional to the parti- 
cle's free Voronoi volume. Thus, the broader the distri- 
bution of Voronoi volumes, the higher the probability of 
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finding larger pores, which leads to a longer tail in the 
exclusion probability as shown in Fig. O The relation 
between the Voronoi volume distribution and the density 
fluctuation follows similar arguments: A broadening in 
the Voronoi volume distribution essentially means that 
there is a broader distribution in the nearest neighbor 
distances and therefore larger differences in the density 
fluctuations. 

We have applied the methods to a set of monodispersed 
spherical particle packings representing stable and coagu- 
lated colloidal particle structures, but the methods could 
of course be generalized. The pore size distribution as 
determined by the Monte Carlo method employed in this 
work can be applied as it is to any porous media. In 
this sense, it is the most general method analyzed in this 
study. The fit using an error function however may not 
necessarily yield good results. In this case the integral 
measure Ips proposed in Sec. IIII Al could be used or the 
width of the distribution could be determined empiri- 
cally. The Voronoi volume distribution generally only 
requires that the elements constituting a structure are 
convex and in this case, the empirical distribution can 
be determined. To the authors' knowledge a fit using the 
/c-gamma distribution has however only been performed 
in the case of packings of monodispersed, spherical par- 
ticles. As for the density fluctuation method, we have 
in this work considered the density of the particle cen- 
ters. The method may be extended to a determination of 
the exact portion of the sphere volumes per cell, which 
however is computationally expensive. An alternative 
could be a cell-wise Monte Carlo integration of the par- 
tial sphere volumes, which would allow for a character- 
ization of arbitrary porous structures using the density 
fluctuation method. 

From an experimental point of view, the methods pre- 
sented in this study rely on the possibility to determine 
the particle positions, which in the case of colloidal parti- 
cles can be obtained using confocal laser microscopy for 
example [24]. In particular, the pore size distributions 
measured using mercury porosimetry [25j and estimated 
using the exclusion probability are not equivalent since 
the latter overestimates the number of small pores due 
to the random probing of the pore space. 

In this paper, we have introduced three scalar mea- 
sures, which for the first time allow quantifying and com- 
paring the heterogeneity of packings of spherical particles 
in terms of a DOH. These measures were calculated using 
distinct techniques and structural characterization meth- 
ods. In view of these differences, the very nice correlation 
between the three DOH-measures is remarkable. Indeed, 
it suggests that the DOH is a microstructure's inherent 
property and that any of the methods proposed in this 
work can be used to uniquely characterize and classify 
it. In terms of sensitivity however, considerable differ- 
ences between the methods where found. Parameter k 
reflects differences in the DOH most sensitively, followed 
by parameter a and finally Idf ■ A further definition of an 
absolute DOH would require a suitable reference struc- 



ture, which for example is either perfectly heterogeneous 
or perfectly homogeneous under the condition of being 
random. 
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